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We study two-particle systems in a model quantum field theory, in which scalar particles and 
spinor particles interact via a mediating scalar field. The Lagrangian of the model is reformulated by 
using covariant Green's functions to solve for the mediating field in terms of the particle fields. This 
results in a Hamiltonian in which the mediating-field propagator appears directly in the interaction 
term. It is shown that exact two-particle eigenstates of the Hamiltonian can be determined. The 
resulting relativistic fermion-boson equation is shown to have Dirac and Klein-Gordon one-particle 
limits. Analytic solutions for the bound state energy spectrum are obtained for the case of massless 
mediating fields. 



I. INTRODUCTION 



o 
o 
o 
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^ ' There are few examples of analytically solvable relativistic two-body wave equations, particularly equations derived 
/••s , from quantum field theory (QFT). The principal exception is the case of spinless (scalar) particles interacting via a 
• mediating scalar field, the so-called scalar Yukawa or Wick-Cutkosky model. This model, therefore, has served as 
a favourite testing ground for methods of solving bound-state problems in QFT, beginning with the Wick Q and 
Cutkosky || solution of the Bethe-Salpeter equation in ladder approximation. 

In earlier papers it was shown that the scalar Yukawa model can be recast in a form such that exact two-body 
eigenstates of the Hamiltonian, in the canonic equal-time formalism, can be determined for the case where there are 
no free (physical) quanta of the mediating "chion" field (i.e. only virtual chions). This is achieved by the partial 
elimination of the mediating chion field by means of Green functions, by the use of the Feshbach-Villars formulation 
■ for the scalar particle fields, and by the use of an "empty" vacuum state. The resulting two-particle bound state mass 
| spectrum, for the case of massless chion exchange, was found to be 



^! E± = \jm 2 +rri$±2m 1 m 2 \/l- — , (1.1) 

where a is the effe ctive dimensionless coupling constant, and n — 1, 2, ... is the principal quantum number. 
q ' The expression (IT) for E(a/n) has th e sh ape of a distorted semicircle. Indeed, in terms of the variables x 



(E 2 — m\ — m|)/2m 1 m 2 and y = a/n, ( |1.1| ) corresponds to half of the circle x 2 + y 2 = 1. The upper branch 



E+ of the energy eigenvalue (1.1) starts at mi + m 2 when a is zero, and decreases monotonically to the value 
E c = y/ m\ + m\ when a reaches the critical value a c = n, beyond which the energy ceases to be real. This 
behaviour is reminiscent of what happens for the one-body Klein-Gordon-Coulomb and Dirac-Coulomb systems. 



X ' I >n r r 2 

From E + = mi + m 2 ^ + 0(a i ) it is clear that E + (a) has the correct non-relativistic Schrodinger 

2 m\ + m 2 vr 

(Balmer) limit. It also has the expected Klein-Gordon one-body limit (with scalar Coulombic potential), namely 
E± — mi — > ±m 2 y/l — a 2 /n 2 as mi — > oo. 

The lower branch, EL, which starts from \m\ — 777.2I, rises monotonically with increasing a to meet the upper 

1 / mim 2 \ a 2 . a 



branch at E c . From £L = Imi — m 2 | H — -, r — r + 0(a ) for mi ^ m 2 , and E— = m h 0(a ) for 

2 \ \mi — m2\J n z n 

777i = = m, it is clear that the lower branch does not have the Balmer limit at small a, and so is 'unphysicaP in 

this sense. The 'unphysical' lower branch arises because an 'empty' vacuum was used, so that positive and negative 

energy solutions (rather than particles and antiparticles) are retained, and this means that solutions with two-body 

energies of the type mi + W2, mi — m 2 , —mi + m 2 , —mi — m 2 occur. However, the use of the empty vacuum is the 

price that needs to be paid for obtai ning analytic solutions. 



Interestingly, the same spectrum (1.1) was obtained previously by Berezin, Itzykson and Zinn- Justin |3| as poles 
of the scattering matrix in eikonal approximation, by Todorov || , using a "quasipotential" one-time reduction of the 
ladder Bethe-Salpeter equation, by Savrin and Troole @ as a solution of the ladder Bethe-Salpeter equation with 
retarding propagators, and recently by Tretyak and Shp ytko Q in a Fokker action formulation. Berezin, Itzykson and 
Zinn- Justin state that their analysis implies that ( |l.l[ ) contains, in perturbation-theory language, all recoil, ladder 
and crossed la dde r effects, but no radiative corrections. However, as has been pointed out recently B], the results 



analogous to (1.1) but for the massive chion exchange case (specifically for fjb/m = 0.15), lie somewhat above the 
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numerical Feynman-Schwinger calculations of Nieuwenhuis and Tjon p0| , which contain all effects save for radiative 
corrections. 

The method of elimination of mediating fields, used in is to some extend similar to that used in the formalism 
of Fokker action integrals in relativistic mechanics. Although the Fokker action formalism deals with a finite nu mbe r 
of degrees of freedom, in the case of scalar interactions its quantized counterpart || leads to the same expression (1.1) 
for the energy as a function of quantum number n as that obtained in . The difference is in different definitions of 
n in the two treatments. This is due to different definitions of scalar interactions in the mechanical and QFT pictures, 
and by quantization ambiguities of the classical Fokker action. 

It is evidently of interest to extend the method employed in to include spinor particle fields. As a first step we 
consider a generalized Yukawa model consisting of fermions, described by a spinor field *, and of bosons, described 
by a scalar field <p, interacting via a mediating scalar field \- O ur starting point is the Lagrangian density (c = fi = 1) 
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C = ^{i^ y d u - gix - mi)* + d v ip*d v p - m 2 ip*ip - g 2 <p*VX + -zpvX&X - 77 A 1 X 



The fields of the model (1.2) satisfy the equations 



(i^dv-giX-mx)^ = 0, 
{d v d v + m\ + g 2X )v = 0, 



and the conjugates of (1.3) and (1.4), as well as 



with 



Equation (1.5) has the formal solution 



(d t/ d v + l x 2 ) X = P, 
P = -(fll** + 52<^» 
X = Xo+Xi, 



Xi = {D*p):= / dx'D(x-x')p(x'), 



where 



dx = d N xdt (in N + 1 dimensions), and Xo( x ) satisfies the homogeneous (or free field) equation, 

(^^ + A1 2 )xo = 0, 

while D(x — x') is a covariant Green function (or chion propagator, in the language of QFT), such that 

[duff* + n 2 )D{x - x 1 ) = 6 N+1 (x - x'). 



(1.2) 



(1.3) 
(1.4) 



(1.5) 
(1.6) 
(1.7) 

(1.8) 

(1.9) 
(1.10) 



Equation (1.10) does not specify D(x — x') uniquely since, for example, any solution of the homogeneous equation 
can be added to it without invalidating (1.1C). This allows for a certain free dom in the choice of D(x), as is discussed 
in standard texts (e.g. refs. jll],[l2]] ) . Substitution of the formal solution (1.7) into Eqs. ( fL3| ) and (1.4) yields the 
"reduced" equations 

{i^d v -g 1 {x +Xx)-m 1 )^ = Q, (1.11) 
{d v d» + ml)y = -g 2 {xv + Xx)v- (1-12) 

These equations are derivable from the action principle 5 / dx C — 0, corresponding to the Lagrangian density 



(1.13) 



C = ^{i'fdu - mi)* + d u ip*d"ip - mlip*<p + xoP + t:p(D * p) 



provided that D(x — x') = D(x' — x). 



The QFTs based on (JE2J) and (1.13) are equivalent in the sense that they lead to identical invariant matri x 
elements in various order of covariant perturbation theory. The difference is that, in the formulation based on (1.13), 
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the interaction term that contains the chion propagator D(x — x') leads to Feynman diagrams that correspond to 
processes involving virtual chions only. On the other hand, the interaction term that contains xo corresponds to 
Feynman diagrams that cannot b e gen erated by the previous term, such as those with external (physical) chion lines. 

The reformulated Lagrangian ( |l.!3| ) contains two types of interactions: "local" interactions, XoP, of the particle 
densities <p*{x)p>{x) and ^>(x)^(x) with the free mediating field Xo(x), and the "nonlocal" interaction, \p(D * p), in 
which the chion propagator appears explicitly. This may seem like a complication rather than simplification of the 
theory based on ( |l.2[ ). However, as we will show, the form (1.13) leads to a model for which exact eigenstates of the 
Hamiltonian can be obtained. 



II. FESHBACH-VILLARS FORMULATION OF THE SCALAR FIELD 



We rewrite the scalar field ip of this model in the Feshbach-Villars (FV) formulation |Tj|] . The reason for doing 
so is that this leads to a QFT Hamiltonian which is Schrodinger-like in form, for which exact eigensolutions can be 
obtained. In the FV formulation, the field ip and its time-derivative ip are replaced by a two-component vector which 
is defined as 



-=i — (m 2 p + ip) 

/2m2 v ' 

(m 2 p - ip) 



/2lTl2 



so that, for example, 2m 2 p>*p> = ((f>l + 2 )(</>i + <fi 2 ) — (ftrrr<f>, where 77 and r are the matrices 



1 

-1 



and t = 



1 1 

-1 -1 



In the FV formulation the equation of motion ( |l.4| ) takes on the form 



i<j) - 



2m 



W 2 T(f) + m 2 r q4> + -J^—Tcfix, 



2m 2 



or, upon using (1.7), the form 



-J—V 2 T(f> + m 2 r ? + ^-T(f>{xo + (D * p)), 
Zm 2 Zm 2 



where p = —(gi^f^f 
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2m,'} 



1 r/rq 



Equations ( |2.4| ) and (1.11) are derivable from the Lagrangian density 
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C FV = ^(i^d v - mi)f + (x)r](j)(x) 
where (j> = ^rjrcj). 



V<f>{x) ■ V0(x) - m 2 </) t (x)(/)(x) + xoP + t:P( d * P)> 
2m 2 2 



(2.1) 



(2.2) 



(2.3) 



(2.4) 



(2.5) 



III. QUANTIZATION 



The momenta corresponding to <p\ and <p 2 are 

BCfv 
P0i = "TT7 = l( Pl 
d(j>i 

Thus, the Hamiltonian density is given by the expression 



(3.1) 



H{x) = ^(x^ix^ix) + ^(x)rjh 2 (x)(t)(x) - xo(x)p(x) - ^p(x){D * p) 



(3.2) 



where hi(x) = —id ■ V + mi/3, h 2 (x) = r ( ) V 2 + m 2 rj, and where we have suppressed terms like 

_ V 2TO 2 y 

V • (cj)(x)V <j){x)) that vanish upon integration and application of Gauss' theorem. We use canonical equal-time 
quantization, whereupon the non-vanishing anticommutation relations, and commutation relations, are 



3 



{V a (x,t),* l3 (y,tf} = 6 a f 3 5 N {x-y), a,0=l..A; [<f> a (x,t),<f>l(y,t)} = r) ab 6 N (x - y), a, 6 =1,2 (3.3) 



where rj a b are elements of the rj matrix (|2.2| ). Using these commutation relations, the Hamiltonian operator can be 
written as 

H = J d N x[H (x)+H x {x)+H I (x)} 7 (3.4) 

where (suppressing the Hamiltonian of the free chion field), 

Hq(x) = * f {x)h 1 {x)^(x) + <f>\x)nh 2 {x)4>{x), (3.5) 
Kx( x ) = ~Xo{x)p{x) (3.6) 



and 



;Hi{x); = -7T / dx' D(x - x') *(x)(*(x / )*(a/))*(a;) 



dx' D(x-x') ~${x) $(x')ct>(x'))<p(x) 
dx' D(x - x')4>(x)(V(x')y(x'))(j>(x) 

dx' D(x - x')^(x)(^{x')(j){x'))^(x), (3.7) 
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.A. 

8m| 

_ gm 

4m 2 

_ gi92 

4m2 

and where we have used the commutation relations ( |3.3|) to reorder 4>{x)4>{x)4>(x')4>(x') as 4>(x)((f>(x')(fi(x'))(f>(x) and 
v E'(x) v E'(x) v I'(a;')5 , (j:') as ^ r (a;)(^ r (a;')^(a;'))^ r (a;). Note that no infinities are dropped upon performing the 'normal 
ordering' of scalar field operators, since none arise on account of the property that r 2 = 0. However, in the case of 
the spinor field, reordering yields an infinite constant which can be absorbed into the total energy of the system. Of 
course, one can simply start from the normal-ordered Hamiltonian. We stress that this 'normal ordering' is not the 
same as the conventional one [T^| , since we normal order the entire field operators, 4>,ip and cft,ip^, not positive and 
negative frequency parts individually. For this reason, we denote it as ; H; rather than as : H :. 

As already mentioned, Tii contains the covariant chion propagator, hence in conventional covariant perturbation 
theory it leads to Fcynman diagrams with internal chion lines. On the other hand, H. x leads to Feynman diagrams 
with external chions. However, we shall not pursue covariant perturbation theory in this work, and so shall not 
consider that app roach further. Rather, we shall consider an approach that leads to some exact eigenstates of the 



Hamiltonian (3.4), but with TL X = 0. 



IV. TRUNCATED MODEL 



In what follows we shall consider the truncated model for which the term 7i x in (3.4) is suppressed. Such a 
Hamiltonian is appropriate for describing systems for which there is no annihilation or decay into chions, or chion- 
phion/psion scattering. 

In the Schrodinger picture we can take t — 0. Therefore, we shall use the notation that, sa y M x, t = 0) = (j>(x), 
etc., for the QFT operators. This allows us to express the interaction part of the Hamiltonian (3.7) as 



■Hi(x); = - Y / d N x'G(x -x')^{x)(^(x')^(x'))^(x) 



8m 2 , 
9i92 
4m 2 

fflff2 

4m? 



d N x' G(x - x')0(x)$(x')(f>{x'))<l)(x) 
d N x' G(x - x') 4>{x)(^{x')t>{x'))<j){x) 
d N x' G{x - x')^(x)($(x')<j)(x'))^(x), 



(4.1) 



where 
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G(x-x') = J^Dix-x'Jdt^j^w I^P^-*'^ 
Explicitly, for N = 3 spatial dimensions this becomes 



M 2 



(4.2) 



for N = 2 it is 



i -»\x-x'\ 

G(x-x') = —- j-, (4.3) 

ill \X - X'\ 



G(x - x') = ^koivlx - x'\), (4.4) 



where ko(z) is the modified Bessel function, whereas for N — 1 it has the form 



G(x - x') = i-e-"'*-*'! (4.5) 



V. EMPTY VACUUM AND ONE-PARTICLE EIGENSTATES 

As in Refs. we define an empty vacuum state, |6), such that 

4>a\0) = * a |0) = 0. (5.1) 

This is different from the "Dirac vacuum" |0) of conventional QFT, which is annihilated by only the positive frequency 
part of <f and \& and by the negative frequency parts of tp' and W . 



With the definition (5.1), the one-particle scalar state defined as 

|1 ) = / d N x<j>Hx)r,f(x)\0), (5.2) 



where f(x) is a two-component vector, is an eigenstate of the truncated QFT Hamiltonian (Ti x = 0) with eigenvalue 
Ei provided that the f(x) is a solution of the equation 

h 2 (x)f(x)=E 1 f(x). (5.3) 

This is just the free-particle Klein-Gordon equation for stationary states (in Feshbach-Villars form). It has, of course, 
all the usual negative-energy "pathologies" of the KG equation. The presence of negative-energy solutions is a 
consequence of the use of the empty vacuum. 
Similarly, the state 

= / d 3 xF a (x)*l(x)\0), (5.4) 



is an eigenstate of ; H; with eigenenergy E\, provided that the four coefficient amplitudes F a (x) are solutions of 

{[hi{x)] a p-E l 5 a0 }F {x)=Q, (5.5) 

or 

(/ii(x) -E 1 )F(x) = (5.6) 
in matrix notation, that is provided the spinor F is a solution of the usual Dirac eigenvalue equation (|5.6|). Note that 



summation on repeated spinor indices a,/3 is implied in equations (5.4) and (5.5) 
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VI. TWO-PARTICLE EIGENSTATES 



We consider a mixed two-particle phion-plus-psion system described by 

|2)= f d N xd N yF aj (x,y)¥ a (x)cj ) \(y)\b) 



(6.1) 

where summation on repeated spi nor and boson indices, a and j, is implied. This state is an exact eigenstate of the 
truncated QFT Hamiltonian ((fo.4|) with Ti x — 0) provided that the 4x2 coefficient matrix F — [F a j] is a solution of 
the two-body equation, 



k(x)F+ \r)h 2 (y)riF T (x,y)\ + V(x - y) 7 ° [r)Tr)F T (x, y)] = EF 

where the superscript T stands for "transpose" . The potential here is given by 

V(x-y) = -^G(x-y) 
2m 2 



(6.2) 



(6.3) 



where G(x — y) is specified in equations (4.2)-(4~5f). Alternatively, we can regard |2) as a variational approximation 
to the eigenstate of the complete Hamiltonian (3.4), since (2\H X \2) = 0. 



For 1 2) to be an eigenstate of the momentum operator with eigenvalue P total = 0, it is necessary that the 4x2 
'hyperspinor' F(x,y) be of the form F(r), where r = x y. Let us define 



where, in N — 3 spatial dimensions, 

s(r) = 



F(r) 



Fu(r) F 12 {r) 
F 2 i(r) F 22 (r) 



s(r) 
t(r) 



t(r) 



F 3 i{r) F 32 {r) 
F 4 i(r) F 42 {r) 



Then equation ( p^2\ ) can be written as two coupled equations for the 2x2-matrices s(r), t(r) 

_ 1 T 



icr ■ \?t(r) + \ r)[h 2 (r) + rV(r) j rj s L (r)\ + (mi - E)s(r) = 0, 
-icr-\/s(r)+ i](ji 2 {r) - rVirfji-jt 1 \r) - (m x + E)t(r) = 0. 



(6.4) 

(6.5) 

(6.6) 
(6.7) 



In the absence of interactions (i.e. V = 0), these equations have solutions with the following four types of energy 
eigenvalues: w(p, mi) +uj(p, m 2 ), —uj(p, mi) +uj(p, m 2 ), to(p, mi) — u>(p, m 2 ), —lu(p, mi)-w(p, m 2 ), where uj(p, m) = 
\l p 1 + m 2 . The first and last are positive and negative energy eigenvalues, respectively, while the middle two are 
'mixed'. The appearance of the negative energies is a reflection of our use of the empty vacuum state, as discussed 
earlier. 



VII. J F EIGENSTATES AND RADIAL REDUCTION IN TV = 3 DIMENSIONS 

The Hamiltonia n (|3.4| ) of the theory commutes with the total angular momentum and parity operators. If the 
two-particle state fl6.1|) is to be an eigenstate of J3 where 



J = J d 3 x^(x,t)j(x)^(x,t)+ J d 3 xcj) i (x,t)r]l(x)(t)(x,t), 



(7.1) 



with j(x) = l(x) + s = —ix x + ^ <T ) then we require that the hyperspinor F (cf. Equation (6T)) must satisfy 
the equation 



ih{x) + l 3 (y))F(x, y) = mjF(x, y) 
In other words, in the rest frame where Ptotal\2) =0 |2), we require that 



(7.2) 
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h(r) + -as j F(r) = m 3 F(r). 
In a similar fashion J 2 \2) = j(j + 1)|2) implies that the matrix F in the rest frame must satisfy the equation 



l' + -+l. (T )F(r)=j(j + l)F(r) 



(7.3) 



(7.4) 



The components s, t of F satisfy Eqs. ( |7.3| ), ( |7.4| ) individually. 

For |2) to be a parity eigenstate, and since V(r) is invariant under space reflection the matrix F(r) must have the 
property that 



j3F(r) = ±F{~r), 
where the ± are the parity quantum numbers. This means that 

s(-r) = ±s(r), t(-r) = Tt{r) 



From Eqs. (|7J) - ([7J) it follows that 
F(r) = 



s(r) 




t{r)_ 





r 

r = -, 

r 



(7.5) 
(7.6) 

(7.7) 



where the upper sign corresponds to parity '+', and the lower to parity and fci, k 2l 9i, 92 are radial functions. 

i—j± - 

The normalized "spinor harmonics" m , 2 (f ) are, explicitly, 



t—i—l I 

C J 2 (f) = 



VW+T 



i + mj + ^Yp 5 (f) 



£_ mi + iy™^(f) 



(7.8) 



and 



C t=3+ Hr)= ' 



VW+T 



(7.9) 



Finally, substituting (7.7) into (6.6), (6.7) we find that the radial functions must satisfy the following system of 
four equations: 

n 2 

9i + 7T^( fc i + fe) + (™i + TO 2)fci + V{r){k! + k 2 ) = Ela, 
2m 2 



j4 

2to9 



(fci + fc 2 ) + (mi - m 2 )k 2 - V(r)(ki + k 2 ) = Ek 2 



n 2 

n a fc i + 77^0?! + 92) + (-mi + m 2 )qi - V(r)( qi + q 2 ) = Eq u 
2m 2 



(7.10) 



n^ 2 



2m- 



■(9l + 92) - (mi + m 2 )9 2 + V(r)(qi + q 2 ) = Eq 2 , 



where V(r) = —a- with a — f 1 ^ 2 ■ 

r %itm 2 ' 



n 



and where we have introduced the operators 



n 2 ± 



d l xl 

r fL + 2 ±- Q' + V2)(j + 1/2±1) 
dr 2 dr r 



(7.11) 
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If we let E — e + mi and consider the mi — > oo limit of (7.10), we find that qi, q% — ► while k\ + k 2 satisfies the radial 
Klein-Gordon equation (with scalar coupling). Similarly, if m\ is replaced by mi in this pro cedu re, then ^2,92 — * 
while q% and — ik\ satisfy the radial Dirac equations (with scalar coupling). Thus, equations ( 7.10] ) have the expected 
one-body limits. 

Defining q± = q\ ± q 2 , and k± = k\ ± fe, then adding the first two, and separately the last two, of equations ( 7.10 ), 
we obtain 



i ~<jr~ = ( E ~ mi ^ k + ~ m2k - + "t 1 T ^ + 5)]?+ 
= (E + mi)q+ - m 2 g_ + ~[1 ± (j + 



(7.12) 



Using definitions (7.11), differentiating (7.12), and substituting the result into ( 7.10| ) we obtain 



2m 2 



(E + mi+ m 2 )q> + m2[lTU + ^ (q ++ q.) 



u- 



- 1 



2r 2 m2 



V 



fc 4 



2rri2r 
m\-\-m-x—E 



(k+- 



2m 2 



2m 2 



2m 2 



(m 2 - mi - E)q' + 



(m 2 — mi + i?)fcl + 



(m 2 + mi - + 



ma[l=F(j + s)] 



2r 2 m2 



2 ' (q+-Q-) 
V 



m 2 [1 ± ( j 



l^2 



1 



2r 2 m2 
m 2 [l±(j + ± 



- F 



(i+|) 2 - 1 



2r 2 m.2 



- V 



k + + 

k-] 
q+ - 
k" 



k' 
2m 2 r 

m\—m<i—E 



(k+-k-) = 



_ 1±Cj+|) , 
2m 2 r ^ + 

-mi+m2-£ 
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2 

Zm 2 r 
mi+m 2 +E 



(q++q-) = 



(?+-?_) = 0, 



(7.13) 



(7.14) 



(7.15) 



(7.16) 



where q' + = etc. Eqs. ( |7.13| ) - ( 7.16 ) do not contain q'_, k'_. Therefore, solving ( 7 . 1 2j ) for g_, putting 

Q + = — ik + , and substituting into (7.13), ( 7.1 5| ) (or into ( 7.141 ), ( 7. 16| ) ) leads to a system of two radial Dirac-like 
equations: 



Q'+ 



-q+ 



i±(i 



m 2 
1^ 



^M<7+ - £ 2<7+ = 0, 



(7.17) 



9+ + "^F(r)Q + - eiQ+ = 0, 



where 



ei = 



(E - mi — m 2 )(E — mi + m 2 ) (£ + mi+ m 2 )(E + mi — m 2 ) 
, £2 = 



(7.18) 



2E 1 2E 

If we make the replacement, E = mi + rn 2 + e, and assume that |V| << mi, m 2 , then equations ( |7.17| ) reduce to 

mim2 



?'j - X 1} Q+ + 2m r ( £ - V(r))Q+ = 0, m r 



mi + m 2 ' 



(7.19) 



which is just the reduced radial Schrodingcr equation for the relative motion of the two particles of masses mi and m 2 , 
with I = j + 1/2 if k = j + 1/2 = £ and I = j - 1/2 if n = -(j + 1/2) = -(£ + 1). Thus, the positive-energy solutions 
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have the expected non-relativistic li mit. On the other hand, for the negative-energy solutions, if E ■■ 
and \e\, \V\ « mi,m 2 , equations ( 7.17 ) reduce to 



q'i- ! ^^q + + 2m r (e + V(r)) q+ 



0. 



-(m 1 +m 2 +e), 



(7.20) 



which is the radial Schrodinger equation, with £ = j + 1/2 if k = —(J + 1/2) = —I and I = j — 1/2 if K = j + 1/2 = £+ 1. 
We note that the sign of V(r) is effe ctively reversed for the negative-energy solutions, relative to that for the positive- 
energy solutions of equation (7.19). This means that if there are positive-energy bound states, then there are no 
bound states for the negative-energy solutions or vice versa. This is similar to that, which occurs for the scalar 
Yukawa model ||,£| . 

We also note that for the mixed-energy solutions of the type E = m\ — m 2 — e, where we take mi > m 2 for 
definiteness, we find that the corresponding non-relativistic reduction leads to an equation exactly like (7.19) but with 
an unphysical reduced mass m r = m\m<i/(m\ — m 2 ). This indicates that (7.17) also admit 'unphysical' mixed-energy 
solutions, in addition to the positive and negative energy solutions. This is consistent with our earlier observation in 
VI regarding the solutions for t he 'free particle' case with V=0. 

The coupled equations (7.17) can be solved readily, for arbitrary mass /x of the mediating scala r fiel d, by standard 
numerical techniques. Once Q+,q+ have been obtained, fc_ and q_ can be determined from (7.12). This, then, 
determines F(r) completely, for any J p state. 



VIII. TWO-BODY BOUND STATES IN 3+1 FOR MASSLESS CHION EXCHANGE. 



We consider the solution of equations ( 7.17 ) for N — 3 and massless chion exchange (i.e. fi — 0), in which 



case the interparticle potential is Coulombic (i.e. V(r) = — ilE 2 . — \ anc i analytic solutions for the eigenvalues and 

87rm2 r 

eigenfunctions can be obtained. Thus, putting 

Q + = e^ r J2 a ^ 1+ ^ q+ = e-f }r J2b j r J - 1+ ' / , (8.1) 

I/=0 !/=0 

we obtain 

(k + 7 + v)a v + ab v — (ia v -\ — £2^-1 =0 (k — 7 — v)b v — aa v + &b v -\ — £\a v -i = 0, (8.2) 



where 



a^, a = -p^-, K = ±(j + \), and 1/ = 0, 1, 2, ... (8.3) 



E 8irm 2 V 2 



The case v — 0, with a u -x = b v -\ = 0, gives 7 = y/n 2 + a 2 . 

In addition, for Q+,q+ to be well behaved at infinity, the series (8.1) must terminate at v = n' > 0. Then (^.2|) , 



with v = n' + 1 and = 6 n '+i = 0, gives the energy spectrum formula 



"'+7=^7^ n' =0,1,2,... where /3 = yf^e^ . (8.4) 
Note that the positive square roots must be chosen for 7 and (3 in order that the wave functions be well behaved at 



the origin and at infi nity. Since n' > and 7, [3 are positive, it follows from (8.4) that a must be positive, which, from 
its definition in (|8.3|), means that for g\g 2 > (i.e. an attractive potential V(r) = —a/r ), the energy eigenvalue E 
must be positive. This, together with the requirement that — t\t 2 > implies that the bound state energy spectrum 
for this fermion-scalar system with scalar Coulombic coupling must lie in the domain \m\ — m 2 \ < E < m\ + m 2 , 



exactly as for the scalar Yuka wa model (1.1) 



It is convenient to rewrite (p.4[) in the form 



ct 2 m% ct , E 2 — m? — mi, 

' + ^VT==^ Wh6re W = 2m im2 - ^ 

The eigenvalue spectrum equation (|S.5|) is not symmetric with respect to the particle masses. This reflects the different 
nature of the particles. The mass mi belongs to the spinor particle and the mass m 2 corresponds to the scalar one. 

Note that by putting n = n' + yjj + 5) 2 + one can rewrite (3.5) in the form 
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E 2 = m\ 



m 2 ± 2m\m 2 



(8.6) 



which formally coincides with the energy spectrum (1.1) for two scalar particles with scalar interaction obtained in 
Refs. [||J3]. But here the 'quantum number' h depends on the energy. 



Equation (8.5) is exact and gives the energy spectrum of the two particle (scalar plus fermion) system. However 
the explicit determination of the function E = E(a,,m\,m 2 ) for various states requires the solution of an algebraic 
equation of the sixth order in E 2 . By contrast, the inverse dependence a — a(E,mi,m2) is relatively simple: 



2m,\E\J\ — w 2 \ 2min'E + yj 'k 2 (E 2 + m\ — m%) 2 + ^m\m\n' 2 (1 — w 2 )^ 



{E 2 



2 

mj 



ml) 



2\2 



(8.7) 



This solution, though analytic, is not particularly transparent, except for some special cases, such as the equal-mass 
states with n' = (i.e. n = n' + j + 1/2 = j + 1/2), for which 



a 



2n\ 



E 
2m 



E = 2m\ 1 



3 + 



1 



(8.8) 



where m — m% = m%. This shows that equal-mass bound states are possible only for a < 2n for n' — (i.e. n — j + |) 
states, and that E = at the critical val ue of a c — 2n. However, for n' > states of the equal-mass case, the shape 
of E(a) is quite unlike the quarter circle fl8.8|) . Rather, E(a) decreases monotonically from E(0) — 2m towards zero, 
as a — * oo. 

As an example, in Figure 1 we plot for n! — states (for which n = j + 1/2), for three different mass 

combinations. The equal- mass curve is labelled mi/7712 = 1, and corresponds to the quarter circle of ( |8.8| ). The 
deformed semi-circle corresponds to m\/m-z — 2. The apex of this curve is the critical value of a beyond which 
there are no real solutions for the two-particle bound state mass E. The two-branch curve straddling the vertical 
asymptote corresponds to mj/mj = 1/2. The vertical asymptote occurs at Efm-i = \J\ — (mi/7712) 2 = V3/2. There 
is a real solution E for any value of a for this mi /mo. < 1 case. Notice that every one of the curves lies in the domain 
I mi — m 2 | <£<mi+ni2 • 
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FIG. 1. Plot of a/n versus Ejmi for n' — (n = j + 1/2) states of three different mass combinations, where mi is the 
fermion mass and 7712 is the boson mass. 



W e can, of course, evaluate and plot a(E) for any n',j state and any values of mi, ma, using the analytic formula 
(3.7). However, it is instructive to outline the general behaviour of E(a), without specifying particular cases. 

First of all, we recall that normalizable bound-state solutions (for which /3 > 0) occur only f or [ mi — m 2 \ < E < 
mi + m.2- There are two branches to the solution E(a), much like for the scalar Yukawa model (1.1) described in the 
introduction. The upper branch, E+(a), begins with the Balmer form E + = m% + m 2 — \ ^'+^2 + 0(a 4 ) at low 
a, then decreases monotonically towards a limiting value E c . There is also a lower 'unphysicaP branch that starts off 



as E_ = |mi-m 2 | + ||^ 
monotonically towards E c 



0(a ) for mi ^ m 2 but as E- = ma/n + 0(a 3 ) for mi = mi = m, and increases 



-ma I 

The lower branch is not of Balmer form at low a. 
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The qualitative behaviour of E(a) is different for the fermion-like case mi < m 2 and the scalar- like case m\ > m 2 . 
This is evident from the one-body limits, described below. For the fermion-like case there are bound state solutions for 
all values of a, no matter how large. Indeed, as a — > oo, the upper branch E+(a) approaches the value E c = yj mi, — m\ 
from above, while the lower branch E-(a) approaches this value from below. Thus, E = E c = yj m\ — m\ is a 
horizontal asymptote of E(a) for the fermion-like cases. 

In contrast, for scalar-like (mi /ma > 1) cases, there are bound state solutions only for finite a < a c , beyond which 
E ceases to be real. The qualitative sh ape of a(E), is that of a distorted upper half-circle (a c being the apex), 
reminiscent of the scalar Yukawa result (1.1). The critical point, E c (a c ) is the end-point for both branches for the 
scalar- like case. The critical value of a varies with mi/ mi > 0. We find that a c /n > 1 for all scalar- like cases. The 
value a c /n = 1 corresponds to the one-body Klein-Gordon limit (mi/m-i — > oo) for which E c — mi = 0. For n' = 
s tate s, a c lies in the domain 1 < a c /n < 2, where a c /n — 1 corresponds to the equal-mass limit, m 2 = mi (cf. Eq. 
(3.8)). For n' > states, a c /n generally increases with increasing n', and becomes arbitrarily large in the equal mass 
limit (mi/m% — > 1). We shall now discuss the limiting cases of the physical, upper branch of E(a) in some detail. 



One mass is large and the other is small 



First, we consider the one-body limits of E + (a), which follow from (3.7) by making the substitutions E — m a (e + l), 

where e = ci-^- + C2 (^fO + and solving for the coefficients c$ . For mi/m-i < 1 we obtain the fermion-like or 
Dirac limit: 



E, 



rn i 

m 2 



< 1 = ni2 + mi 



'1- 



mi 



m 2 2(ri + Va 2 + k 2 ) 2 



(vf + Va 2 + k 2 ) 2 

2o^ 

yju 2 + n 2 (ri + Va 2 + k 2 ) 




(8.9) 



The first term in the square brackets is just the Dirac one body energy spectrum for a Coulombic potential (scalar 
coupling), and the second one gives the first order correction in mi/m.2. 
For mi I mi < 1 the expansion yields the scalar or Klein-Gordon limit: 



^4 



m 2 
mi 



< 1 = mi 



m 2 




(8.10) 



The first term in square brackets coincides with the Klein-Gordon one body energy spectrum in a Coulombic potential 
(scalar coupling), while the second gives the first order correction to it, in powers of m 2 /mi. 

Expansion of E in powers of the coupling constant 



Next, we consider the expansion of E + (a) in powers of a. The result is 



1 



1 



E + (a>) = mi + m 2 m r — =■ m r a 



mim 2 



-i-4- 



+ 0(a 6 ), (8.11) 



(mi+m 2 ) 2 /n 4 (mi + m 2 ) 2 n 3 (j + |) _ 
We obtain the expected non-relativistic Balmer result at 0(a 2 ). The 



, 1 mim 2 

where n = n + j + — and m r = . 

2 mi + m 2 

0(a 4 ) correction is not symm etric in mi and m 2 due to the different, fermionic and bosonic, nature of the particles. 



The one body limits of (3.11) have the required Dirac and Klein-Gordon forms, as can be seen from the expressions 



£4 



mi 
m 2 



< 1 = m 2 + mi 



1 



1 



+ 



mi / a 2 

m~ 2 ~ 2n 3 (j + i) 



n 4 n 3 (j + l) / 
3a 4 x / /_ x 2 



i? + ( — < 1 
mi 



mi + m 2 



In 2 8n 4 mi 2n 



Note that (3.12) and (p. 13) are the same as the expansions of Eqs. (|8.9|), (p.lO) in powers of a, as they ought to be 




m 2 a 



.12) 



(8.13) 
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IX. CONCLUDING REMARKS 



We have studied two-particle systems in a model QFT, in which fermions of mass mi interact with bosons of mass 
mi- The interaction is mediated by a real, scalar field of mass fi (the 'chion' field). The field equations were used to 
recast the Hamiltonian of the theory into a form in which the chion propagator appears directly in the interaction 
term. For the case where there is no decay, emission or absorption of real (physical) chions (i.e. only 'virtual' chions), 
we obtain exact two-particle eigenstates of the Hamiltonian, using an unconventional 'empty' vacuum state, which is 
annihilated by both the positive and negative frequency parts of the particle-field operators. 

The resulting relativistic two-particle wave equation, for the stationary states of the system, reduces to a pair 
of Dirac-like radial equations for the various J p states. These equations are shown to have the radial Schrodinger 
equation for the relative motion of the two particles as the non-relativistic limit, and the Dirac and Klein-Gordon 
equations (with scalar coupling) as the one-body limits. Analytic solutions for the two-body bound state eigenenergies 
(rest masses) are obtained for the massless chion exchange (fi = 0) case. The shape of the E(a), or a(E), curves, 
where a is the dimensionless coupling constant, is discussed for various mass combinations, mi/m2, and various nJ p 
states. 

In the case of massive chion exchange (fi ^ 0), the eigenvalues and eigenfunctions must be obtained numerically, 
which can be done easily by standard methods. We do not present such solutions in this paper. Also, we do not 
discuss the scattering-state solutions of the equations, though these can be worked out readily. 

Lastly, we mention that A/"-body eigenstates, where Af > 3, and the corresponding relativistic A/"-body equations, 
can be worked out readily for the present model, as was shown for the purely scalar model in || , and for QED in [|l5| . 
Such equations are much more complicated, since they possess all the complexity of relativistic many-body equations. 
We do not discuss the J\f > 3 systems in this paper. 
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